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ABSTRACT 



The extreme, intra-hour and > 10% rms flux density scintillation observed in AGNs such as PKS 0405-385, J18 19+3845 and PKS 
1257-326 at cm wavelengths has been attributed to scattering in highly turbulent, nearby regions in the interstellar medium. Such 
behavior has been found to be rare. We searched for rapid scintillators among 128 flat spectrum AGNs and analyzed their properties 
to determine the origin of such rapid and large amplitude radio scintillation. The sources were observed at the VLA at 4.9 and 8.4 
GHz simultaneously at two hour intervals over 11 days. We detected six rapid scintillators with characteristic time-scales of < 2 
hours, none of which have rms variations > 10%. We found strong lines of evidence linking rapid scintillation to the presence of 
nearby scattering regions, estimated to be < 12 pc away for ~ 200 //as sources and < 250 pc away for ~ 10 //as sources. We attribute 
the scarcity of rapid and large amplitude scintillators to the requirement of additional constraints, including large source compact 
fractions. Jl 8 19+3845 was found to display ~ 2% rms variations at ~ 6 hour time-scales superposed on longer > 11 day variations, 
suggesting that the highly turbulent cloud responsible for its extreme scintillation has moved away, with its scintillation now caused 
by a more distant screen ^ 50 to 150 pc away. 
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1. Introduction 

Hourly variations in the flux densities of the most compact active 
galactic nuclei (AGNs) at cm wavelengths have incontrovertibly 
been shown to result from interstellar scintillation (ISS). The ev- 
idence comes from observations of time delays between vari- 
ability patterns across widely- separated telescopes ([Jauncey et 



al. 20001 Dennett-Thorpe & de Bmyn 2002| BignaU et al. 2006), 



as well as annual cycles in their variability timescales (Rickett 
et al. 2001; Jauncey & Macquart 2001 [ [BignaU et al 



only a small fraction of the sky. It also begs the question as to 
whether the presence of these nearby clouds are the only neces- 
sary condition for rapid, large amplitude ISS, or if there are other 
factors that contribute to their scarcity (ie. source compactness). 



Jl 8 19+3845 was the most extreme of this class. During a 
4.8 GHz, 6-hour monitoring program on the European VLBI 
Network (EVN) in 2008 March, its extreme scintillation sus- 



^QQ^I tained over at least eight years was found to have ceased ([Cimo 



Dennett-Thorpe & de Bruyn,20Q3, Jauncey et al.|2QQ3| ) consis- [p08 ), with an upper limit of 1% rms variations (Cimo, private 
. . communication). Understanding the unexpected disappearance 

of the extreme variability of this source will shed some light on 
its origin at earlier epochs. We note that epis odic extreme scintil- 
lation has been observed in PKS 0405-385 ( |Kedziora-Chudczer| 
|2006 ), attributed either to the repeated appearance and expan- 
sion of new components in the source or spatially intermittent 
turbulent patches in the ISM drifting across its sight-line. 



tent with relative motion between the scattering region and the 
Earth. These observations were made possible through the mon- 
itoring of the so called 'extreme scintillators' that display rms 
flux density variations > 10% on timescales < 1 hour, of which 
PKS 0405-385 ([Kedziora-C hudc zer et alTJTWf) , J1819+3845 
([Dennett-Thorpe & de Bru yn 2000 ) and PKS 1257-326 ( [BignaU 



et al.|2Q03| ) are the most well-known. Such rapid and large am- 



plitude ISS has been attributed to the presence of nearby turbu- 
lent clouds, < 30 pc from the Sun. 

These extreme scintillators are surprisingly rare, considering 
that the three archetypal sources were detected serendipitously. 
In the 5 GHz, Microarcsecond Scintillation Induced Variability 
(MASIV) Survey of 443 compact flat- spectrum sources ( Lovell, 
et al. [[2008] ), no new examples (apart from Jl 8 19+3 845) were 
found to display sustained modulation indices (rms variations as 
a percentage of mean flux density) of more than 10%, although 
16% of the sources were found to scintillate on timescales of 
< 12 hours. This implied that the nearby, turbulent clouds cover 



In this letter, we examine data from a dual-frequency survey 
of ISS in 128 flat- spectrum AGNs to search for rapid scintilla- 
tors, then make use of their properties to understand the origin 
of rapid radio scintillation. We present a method for identifying 
rapid scintillators and include a list of candidates in Section [2| 
In Section [3] we report on the variability of Jl 8 19+3845 sub- 
sequent to the period of extreme scintillation. We then discuss 
the implications of our results on the physics behind the most 
extreme scintillators in Section (4] 
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Fig. 1. Do/,^(2day) - Dnoise plotted against Z)o^^(2hr) - Dnoise at 
8.4 GHz (top) and 4.9 GHz (bottom). For each source, only the 
lower limits of the Icr errors in D{2 day) and the upper limits of 
the Icr errors in Z)(2hr) are shown. The solid lines have a slope 
of 24, 2 and 1 (from top to bottom) respectively. The dashed 
lines represent the 5cr thresholds. Rapid scintillators are located 
within the shaded regions. 



2. Detection of the Six Rapid Scintillators 

As a follow-up to the MASIV Survey, a sub-sample containing 
128 of the 443 origina l MASIV sources was observed at the 

hereafter Paper I). 



VLA in 2009 January (Koay et al. 2011 



The increased sensitivity from the larger number of antennas per 
subarray and observation span of 1 1 days enabled the measure- 
ment of rms flux density variations down to a level of ~1% and 
the estimation of scintillation timescales. Simultaneous obser- 
vations at 4.9 GHz and 8.4 GHz on two subarray s provided a 
robust means of discriminating flux variations of astrophysical 
origin from variations due to instrumental and systematic errors. 

We characterized the variability of each source using the dis- 
crete structure function (SF) of the lightcurves, given by: 



Dobsir) 



J,k 



where S(t) is the flux density of the source at time t, normal- 
ized by its mean flux density. A^^- is the number of pairs of flux 
densities with a time-lag r. We also estimated a constant Dnoise 
for each source, which accounts for variability arising from in- 
strumental and systematic errors including noise, low-level con- 
fusion, as well as residual gain and calibration errors, which we 
then subtract from Dobsij) to obtain the 'true variability' of each 
source, Dij). The details are described in Paper I. 

We found significant variability in 30% of the sources at 
2-hour time- scales, the shortest time interval between succes- 
sive observations of each source. These sources have D{2 hr) > 
2 X 10"^ on at least one frequency, so that DobsC^ hr) > 3cr above 
Dnoise (~ 1 X 10""^ on average at both frequencies). These sources 
also display strongly correlated variability at both frequencies, 
based on the cross-correlation function, and a visual inspection 
of the lightcurves. Such correlations are unlikely to result from 
antenna-based or array-based errors as they were observed on 
two separate subarray s. Neither are these variations likely to be 
confusion, since the subarrays each have a diff'erent uv-coverage 
arising from dissimilar primary beams and synthesized beams at 
the two frequencies. On the other hand, flux variations due to 
weak ISS are expected to be correlated across a wide-bandwidth 
(NarayanJ992) . We also found larger fractions of sources dis- 
playing significant 2-hour variability at 4.9 GHz and in the weak 
sample of sources (Sy < 0.3 Jy), consistent with a population of 
brightness temperatu re limited sources w hose variability is dom- 
inated by weak ISS ( Lovell et al. ||2008| ). No similar variations 
were observed in sources close to each other, ruling out atmo- 
spheric eff'ects. None of the sources were observed at low solar 
elongations, so interplanetary scintillation is also negligible. 

We devised a method to distinguish the rapid scintillators 
with characteristic timescales of Tchar < 2 hours from the sources 
with much longer Tchar that display significant 2-hour variability. 
We have defined Tchar as the time-scale at which D(t) reaches 
half its value at saturation. D(t) increases with and saturates at 
twice the true variance of the source, when the source is observed 
for a duration much longer than Tchar- Letting Tsat be the time- 
scale at which D(t) saturates, we obtain: 



if Tsat < 2 hr 

if 2 hr < Tsat < 2 days 

if Tsat > 2 days 



(2) 



S(tk-T)r 



(1) 



^^| = (W2hr)« 
D(2hv) |>24f 

Figure [T] plots Dobsi^day) - Dnoise against Z)o^^(2hr) - Dnoise 
at both frequencies so that Z)(2day)/Z)(2hr) is represented by 
a slope. We determined that a source must satisfy both of the 
following conditions at at least one frequency to be classified as 
a rapid scintillator, corresponding to shaded regions in Figure [T] 

1. Z)(2day)/Z)(2hr) < 2, since Tchar > 0.5(t,«^) for a > 1. 
We include also sources that lie outside this region, but have Icr 
error bars that extend into this region. 

2. D(2 day) and D(2 hr) must be > 4 x 10""^, so that Dobs(r) > 

5cr above Dnoise- 

We detected six rapid scintillators in the sample (listed in 
Table [T]), out of an initial ten candidates that fulfilled the con- 
ditions above. Four candidates were found to be slower scin- 
tillators after a visual inspection; they contain outliers in the 
lightcurves that skew D(2hr) towards larger values, or dis- 
play quasi-periodic structure in the lightcurves leading to quasi- 
periodicity in the SFs, thereby reducing D(2 day) when 2 days is 
close to a multiple of the quasi-period. For the six confirmed 
rapid scintillators, a cross-check with their Tchar obtained us- 
ing an independent method in Paper I revealed that they were 
also estimated to be < 2 hours on at least one frequency. 
The lightcurves of these six sources show strong correlation 
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SOURCE: J1 328+6221 



SOURCE: J1 81 9+3845 
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Fig. 2. Lightcurves for the rapid scintillator J1328+6221 at 8.4 
GHz (top) and 4.9 GHz (bottom). 
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Fig. 3. Lightcurves for the former extreme scintillator 
J1819+3845 at 8.4 GHz (top) and 4.9 GHz (bottom). 



Table 1. Rapid Scintillators Detected in our Sample and Their 
Properties. 



Source Name 


Galactic Coordinates 


4" 


o b 

^4.9 


c b 

^ 8.4 






(J2000) 


/ 


b 


(R) 


(Jy) 


(Jy) 


(%) 


(%) 


J0800+4854 


170.11 


31.16 


0.5 


0.10 


0.08 


5.0 


5.1 


J0929+5013 


167.66 


45.42 


0.6 


0.40 


0.39 


4.1 


2.6 


J1049+1429 


230.79 


59.00 


0.9 


0.13 


0.15 


3.0 


2.8 


J1328+6221 


115.56 


54.23 


0.5 


0.10 


0.08 


5.9 


3.0 


J1549+5038 


80.24 


49.06 


0.4 


0.91 


0.93 


3.0 


2.1 


J 193 1+4743 


79.98 


13.53 


5.2 


0.11 


0.10 


7.2 


7.4 



Notes. Line-of- sight Ha intensity in units of Rayleighs (R), from 



Haffneretal!l ( [2003 ) 2009 January VLA mean flux density 2009 



January raw modulation index, without error subtraction 



across both frequencies. An example, J 1328+6221, is shown in 
Figure [2] 

3. The Variability of J1 81 9+3845 

Our observations confirm that the character of the rapid and large 
amplitude scintillation in Jl 8 19+3 845 has changed significantly 
(it is no longer exhibiting extreme scintillation, see Figures[T]and 
|3]), as previously noted by Cimo_ ( 2008 ). We detected rms varia- 
tions of only 1.5% at 4.9 GHz and 2.5% at 8.4 GHz. Variations 
of Tchar ~ 6 hours dominate at 4.9 GHz, while the variability at 
8.4 GHz is dominated by the slow rise in flux density over the 
11-day period, also discernible in the 4.9 GHz lightcurve. Such 
variations would not have been detected in the 6-hour observa- 
tions by Cimo. 

There are two possible explanations for the factor of ~ 20 de- 
crease in scintillation amplitude and the factor of ~ 12 increase 
in timescale from the previously observed rms variations of 25- 
40% and time-scale of 30 minutes ( Macquart & de Bruyn 2007 ). 
Either the source has expanded, thereby quenching the scintil- 
lation, or the highly turbulent patch in the Local Interstellar 
Medium (LISM) responsible for the extreme scintillation has 
drifted ofl' its sight-line. If the former is true, it would require the 
apparent angular size of the source to have expanded by a factor 
of 12 in < 2 years from 2006 to 2008. For a > lOOyuas source 
(65 pc at its measured redshift of 0.54, Hq = 70kms"^Mpc"\ 
flM = 0.27 and ^Ia = 0.73), this requires an apparent expansion 
speed of > 5.8c. As there is no discernible change in its mean 
flux density, and the source remains unresolved at mas sc ales at 
all EVN baselines which includes the Urumqi telescope ([Cimo 
2008| ), such an explanation is unlikely. 



The latter explanation is therefore more likely. Assuming 
that the source size has not changed and that the velocity of the 
screen is 50kms"\ we estimate that the scattering region now 
lies at a distance of 50 to 150 pc away, requiring also a factor 



~ 9 decrease in the length-scale, rdiff , of the phase fluctuations 
at the scattering screen to achieve the observed reduction in the 
modulation index (see Equation |3]). The slower, Tsat > H day 
variations could be intrinsic to the source, or IS S at a second, 
more distant screen -1.7 kpc away. In any case, these longer 
time- scale variations would have been masked by the previous 
rapid and large amplitude variations of Jl 8 19+3845. 

4. The Origin of Rapid and Extreme Scintillation 

The spatial distribution of the six rapid scintillators shows a de- 
pendence on Galactic latitude and line-of- sight free electron con- 
tent, consistent with the view that the most rapid ISS in quasars 
is caused by scattering in nearby regions. Of the six sources, five 
of them lie at Galactic latitudes of > 30°, and have Ha intensi- 
ties obtained from the Wisconsin H- Alpha Mapper Northern Sky 
Survey ( [Haff'ner et al.|2003| of /q, < 1.0 Rayleighs, as a proxy for 
the line-of- sight integral of the square of the free electron den- 
sity. In fact, out of 20 sources in our sample that have la < 0.6 
Rayleighs, of which 14 have variability amplitudes 3cr above 
D noises four of them are rapid scintillators. On the other hand, 
in 108 sources with I a > 0.6 Rayleighs, we find only two rapid 
scintillators. Fisher's Exact Test for contingency tables confirms 
that the association between rapid ISS and low la is statisti- 
cally significant at the 0.05 level. We calculated a probability of 
0.0044 that the correlation was obtained purely by chance when 
considering the entire sample of 128 sources, and a probability 
of 0.0026 when considering only the > 3cr variable sources. The 
significance holds even if we consider Jl 8 19+3845 (7^, = 2.2 
Rayleighs) as a rapid scintillator based on its history, and include 
PKS 0405-385 and PKS 1257-326 in our sample, which have 
/q, < 0.5 and 7^, = 16 Rayleighs respectively (obtained from the 
Southern H- Alpha Sky Survey Atlas, see Gaustad etal. "(2001)). 
This is consistent with the MASIV results which show that the 
fraction of fast scintillators ijchar < 0-5 days) decreases with 
increasing Ha intensity ("L ovell et al. |2008| ). These sight-lines 
with low Ha at higher Galactic latitudes correspond to regions 
of lower transition frequencies (vq) between weak and strong ISS 
( Walker,^ 1998) , which imply lower eff'ective scattering screen 
distances at a fixed observing frequency. These nearby screens 
produce smaller scale interference patterns on the Earth's sur- 
face, leading to more rapid scintillation for a given screen veloc- 
ity as they drift across the telescope. 

The Tchar < 2 hours can be attributed to scattering at an ef- 
fective screen distance of L < 12 pc for a ~ 200yuas source com- 
ponent, or at L < 250 pc for a ~ lOyuas component, for typical 
transverse velocities of 50kms"^ Assuming Kolmogorov turbu- 
lence and vo = 3 G Hz, numerical computat i ons us ing the fitting 
formula for ISS in Goodman & Narayan ( 2006| ) indicate that 
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the former will produce the observed ~ 5% rms variations, giv- 
ing brightness temperatures, Tg, of 4 x 10^^ K if all the 100 mJy 
flux density is contained within the ~ 200 yuas component. The 
latter requires that the ratio of the ~ lOyuas component flux den- 
sity to the total observed flux density (ie. the source compact 
fraction, fc) be ~ 25% to obtain -5% rms variations, giving 
Tb = 4 X 10^^ K for an observed 100 mJy source. This is be- 
cause the observed modulation index is normalized by the flux 
density measured at the resolution of the telescope, which may 
include emission from extended, non- scintillating components. 

There are strong indications that both scenarios for rapid ISS 
occur in our sample. Gaussian mo del fits to so urce core com- 
ponents in 8.4 GHz VLBI images ( Qjha et al. 2004) estimate 
core diameters of ~ 200 //as (L < 12 pc) in J1328-F6221 and 
J1931-h4743, as weU as ~ 40yuas (L < 60 pc) in J1049-hl429. 
The VLBI cores may contain an even more compact, unresolved 
component which could push the maximum screen distances fur- 
ther out. However, |Linsky et al. ( 2008) have also shown that the 
extreme scintiflation in J1819-F3845, PKS 0405-385 and PKS 
1257-326 may be associated with regions where multiple warm- 
ionized clouds < 15 pc from the Sun intersect, and possibly inter- 
act to form highly-turbulent screens. We found that J1328-h6221 
and J1931-h4743 also have sight-lines through (or at the edge 
of) these multiple clouds, providing further evidence for screens 
< 15 pc. On the other hand, the sight-lines of J1049-hl429 and 
the other three sources do not intersect these regions, suggesting 
larger screen distances. 

That rapid ISS arises from a large range of source sizes and 
scattering screen distances explains why rapid scintillators are 
not as rare as initially thought; it is the sources that display both 
rapid and large rms variations > 10% that are rare. The six rapid 
scintillators constitute 5% of our sources, yet none have rms 
variations > 10%. No new extreme scintillators were found in 
the MASIV Survey ( jLovefl et al. f 2008 ), although J0949-h5819 
and J0829-h4018 did show ~ 15% rms variations at one epoch. 

We argue that the scarcity of the extreme scintillators can be 
attributed to additional constraints required for rapid scintilla- 
tors to also display rms variations > 10%. In the regime of weak 
ISS where the observing frequency, v > vq, the modulation in- 
dex for an extended source whose angular size, Og, is larger than 
the Fresnel ang ular- scale at the sca ttering screen. Of, is given by 
Narayan (1992]) and Walker (1998]), which we rewrite as follows: 



27ry\rdiff 



Os 



(3) 



where Of ^ r^lL = ^/c/(2nyL), rp is the Fresnel length-scale, 
rdiff is the difl'ractive length-scale determined by the strength of 
the turbulence at the scattering screen, and rp/ rdiff = (vo/v)^^^^^. 
Kolmogorov turbulence is assumed. 

Firstly, the source must be observed at a frequency close to 
vo, where m for a point source -100%. This is also true in the 
regime of strong refractive ISS (v < vq), as m decreases again 
with decreasing y. 

Secondly, the scattering screen must be highly turbulent 
(with small rdiff). As discussed in Section |3] such a turbulent 
cloud has possibly moved away in the case of J1819-h3845. 
These turbulent patches may be localized and intermittent, pos- 
sibly causing the episodic extreme scintillation in PKS0405-385 
( [Kedziora-ChudczeF^ |2006| ). Such AU-scale inhomogeneities 
deviate from standard models of the ISM with homogenous 
Kolmogorov turbulence. They are instead reminiscent of the 
clumpy, discrete clou ds observed in extreme scattering events 
( [Fiedler et a'T~||1987| ) and the scintillation of pulsar B0834-h06 



(| Brisken et al. |2010] ). Scattering in such clouds would dominate 
the scintillation in PKS 1257-326 over any background ISS, de- 
spite its large line-of-sight 7^,. 

Thirdly, the source must have fc close to unity. Comparing 
VLBI flux densities of the core co mponents to that of the ex- 
tended mas structures (Qjh a et aLl [2004) gives upper limits to 
fc at VLBI scales; < 75% for J0800-h4854, J1328-h6221 and 
J1931-h4743, < 100% for J 1049 -hi 429. J1549-h5038 also shows 



significant extended structure at mas scales ( |Fey & Chariot 
(2000). At the resolution of the VLA, fc will possibly be lower. 
This reduces th e modulation indices in our six rapid scintillators. 
As an example, Savolainen & Kovalev (2008 ) reported m ~ 13% 
scintillation at a time-scale of 2.7 hours in J1159-h2914 with a 
mean flux density ~ 1.5 Jy at 15 GHz, based on VLBI observa- 
tions in 2007. The MASIV Survey and our observations found 
m ~ 5% at both frequencies, obtaining flux densities ~ 3.0 Jy. 
The unnormalized peak to trough variations in the three separate 
experiments are in fact comparable in amplitude, ^ 0.5 Jy. 

In summary, while scattering screens at distances up to 250 
pc can produce Tchar ^ 2 hour ISS in AGN cores > 10 yuas, source 
compact fractions > 50% are necessary for ~ 10 yuas scintil- 
lating components (> 95% for -100 yuas components) to dis- 
play > 20% rms variations at Tchar ~ 1 hour when observed at 
y = Vq ~ 4.9 GHz. We therefore predict that higher angular reso- 
lution observations will reveal more large amplitude rapid scin- 
tillators, providing more reliable constraints on the flux density 
of the scintillating components. 
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